THE MODEL MATCHING PROBLEM FOR A GENERAL CLASS OF NON- 
LINEAR DISCRETE INPUT - OUTPUT SYSTEMS WITH CROSS-PRODUCTS. 
AN ALGORITHMIC APPROACH. 
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Abstract: In this paper we examine the model matching problem that concerns nonlinear input - 
output discrete systems, containing products among delays of input and output signals, through a 
special factorization. The algebraic framework of (5e-operators and the star-product that we adopt 
describe these systems. Moreover, a certain procedure that resembles the Euclidean division, 
\ allows us to discover the linear factors of those systems, with respect to the above mentioned 

^ ' operations. The entire approach is symbolically algorithmic and involves the use of suitable 

software. 
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As is well known, nonlinear systems are used in a variety of applications and have been 
the focus of research for a number of years. This, along with a recent increase in their use, 
has led to a steady development of the relevant theory. Furthermore, various advances in 
computer technology have had a major effect on control system analysis and design. These 
advances have led to the development of new algorithms and techniques, which in turn 
have allowed design methods to be accomplished with greater speed and efficiency, [13]. 
^ I Yet, computational methods for the study of nonlinear systems are (still) at an early stage 

^ ■ of development. Such computational methods have been applied in the study of modeling, 

the problem of feedback linearization, and in the global optimization problem, ([3], [4], [13], 
to mention but a few). All of these problems were investigated mostly in relation to 
continuous systems, and, in some cases, in relation to non-linear discrete time systems, [14] . 
By definition, a discrete time system evaluates input and output signals over a countable 
number of time instants. In some cases, a discrete time system is obtained from continuous 
time systems through sampling at certain time instants [17]. A further case concerns 
systems that are naturally and directly described in discrete form, typically in financial 
or economic systems. There is a rich literature devoted to the study of discrete systems. 
Certain works approach the issue through analytical tools, as the Lie Derivatives and 
Distributions [6], [8] and others by using algebraic methodologies, like differential algebra 
or rings theory, [3], [4], [5], [7], [27]. Some researchers give algorithmic results [5] and others 
study the discrete systems from the theoretical point of view [23], [24] An important feature 
of discrete systems is their strong computational orientation. This is connected to the fact 
that the procedures developed for these systems are more " suitable " for realization in 
computational machines. Thus, several design techniques can be accomplished at a faster 
rate and in a way that makes them appropriate for further applications. 



The present paper focuses on nonlinear discrete input-output systems of the form: 

y{t) + Yl ^ivi* -^) + J2Yl (^hnvi^ - '^i)y{t - ^2) h \- 

= Yl - j) + --- + Y12---Yl ^hh-3mU{t - Ji) ■■■u{t- + 

+EE^^^?/(*-^M^-0 + --- + EEE---EE---E 

Ckik2...K,hi2...ivy{t -ki)--- y{t - K')u{t -h)--- u{t - Q (1) 

Equation (1) transforms causal input signals (i.e. u{t) = for i < 0) to causal output 
signals. A set of initial conditions j/(0) = |/o, yil) ~ yi, ■ ■ ■, y{k) = Uk is always assigned to 
(1) and the lowest delay output term (that is y{t)) appears in the linear part of the system. 
From all the above, wc conclude that for any given input sequence the entire future output 
is uniquely determined, [16]. The systems of the form (1), which contain products among 
inputs and output signals, that are sometimes called 'cross-products', encompass a broad 
variety of nonlinear discrete systems. Wc obtain these either through transformations of 
nonlinear discrete state-space representations into input - output forms [20], or when we use 
Taylor's expansion method to approximate other more general nonlinear discrete systems 
[16], [20]. These are employed in signal processing theory, whenever it is necessary to 
construct nonlinear representations of discrete signals, (they are an extension of the infinite 
inpulse response filters to a nonlinear set up, [16]), in nonlinear time-series analysis and in 
adaptive control, in the context of designing nonlinear adaptive controllers, [26]. Design 
problems for such systems have been studied in the past through the use of cither analytical 
or algebraic tools (see for instance [3], [5], [19]). We have to mention here especially the 
following works [23], [24] which study the above systems by means of algebraic methods 
and provide us with theoretical results. 

In this paper wc would like to propose a novel way of dealing with systems of the form (1) 
within a computer algebra environment. We do this initially by means of special operators 
and operations that create an appropriate algebraic framework and then through the devel- 
opment of symbolic computational algorithms that allow us to solve the model matching's 
design problem. Specifically, we use the so-called 5e-operator to deal with cross-products. 
This operator, introduced in [10], is an extension of the simple S and e-opcrators that 
arc used in the algebraic description of nonlinear discrete input - output systems without 
cross-products [11], and which are also an extension of the simple shift operator q in linear 
discrete systems [21]. Via these operators we can define the so-called 5,e and Se - poly- 
nomilas. These polynomials involve the delays of output and input signals in a nonlinear 
polynomial way, and are similar to polynomials with many variables. ( A special case of 
them are the linear polynomials ). By means of these polynomials we can formally rewrite 
(1) as follows: F[y{t),u{t)] = 0, where F a 5e-polynomial. Among these polynomials, we 
can define two product operations: the dot-product, which corresponds to the usual prod- 
uct among polynomials, and the star-product, which corresponds to the substitution of one 
polynomial by another, or, in more system-oriented terminology, to the cascade connection 
of systems. Though there is a similarity between this algebraic background and others, 
there is also the interesting peculiarity that the set of 5e-polynomials with respect to the 
star- product is not a ring, [3]. Indeed, the property A* (B + C) — A*B + A*C is NOT 



The concept of formal 5 or e-linear factors is presented here and departs from previous 
approaches, [9]. The basic idea underlying this concept, is that it " reveals ", through 
appropriate algorithms, the linear S or e-polynomials that are hidden within a nonlinear 
(5e-polynomial F and they are factors of F with respect either to the dot-product or to the 
star-product. Specifically, if F is a given 5e-polynomial, then we can rewrite it as: 

F = Y: S^o'^o' ■ {ceSi^e^e * [Le, Me]) + Rs,h + "^e^ + %e (2) 
e 

where, in general, the coefficients cq, the linear polynomials L^, Mq and the remainders 

^,Re_/i,R^g, depend on the parameters Wijhk,Sijhk- These polynomials Lg,Mg arc the 
linear formal factors of F, upon discussion. The immediate consequence of this factorization 
is that we can transfer questions about the behaviour of the nonlinear system to questions 
about the structure of the linear polynomials it " contains " and apply well-known methods 
for dealing with those polynomials. The above expression of a given 5e-polynomial, is 
achieved by means of certain computational algorithms. These algorithms are simply a kind 
of division among 5e-polynomials with respect to the star-product. Though this procedure 
resembles the method of Grobncr basis [2], there arc some essential differences: (1) The set 
of 5e-polynomials is not a ring. (2) Here we are using just one 5e-polynomial, instead of a 
set of such polynomials. (3) It involves the notion of free parameters. We can then achieve 
certain tasks by giving suitable values to those parameters. For instance, we can eliminate 
the remainders R^ ^, Rg./j, R^^ or we can find the common factors of the polynomials 
LqjMq, if any. All these can be implemented in a computer and an appropriate software 
has already been created through the Mathematica Package, which was used throughout 
this paper. To elucidate the whole issue, let us present a short exhibitory example. Let 
us have the bilinear system: y(t) + y{t — 1) = u{t — 1) -|- y{t — 2)u(t — 1). This system 
can rewritten as {Sq -t- — ei — 62ei)[y{t),u{t)] = or F[y{t),u{t)] — 0, where F is the 
5-polynomial Sq + Si — ei — S2ei. F can be written as: 

5o + 5i - €i - 5oeo * [wi,o,o,o(^o + «^i,o,o,i^i + (^2, Si,o,o,oeo + ei] + Si,o,o,o^oeo * [«^2,o,o,o^o+ 

+«^2,o,o,i^i + ^2, eo] + «^i,o,o,i<^oeo * ['«^3,o,o,o'^o + ^1, 53,0,0,0^0 + ei] + («^i,o,o,iSi,o,o,o- 
-«^2,o,o,iSi,o,o,o - '"^1,0,0,1 S3,o,o,o)^oeo * [w^4,o,o,o^o + ^1, eo]+ 
+(«^i,o,o,o - «^i,o,o,i^y3,o,o,o)(^oeo * [So, 55,0,0,0^0 + ei] + R^^ 

(we do not write R^^ explicitly, due to its large size). The exact meaning of the indexes 

in the above expression will be explained later. The set of the Formal Linear factors of 
F is: {60 + Si,ei, Wi,o,o,o<^o + «'i,o,o,i'^i + (^2, Si,o,o,oeo + Ci, ^2,0,0,0(^0 + ^^^2,0,0,1(^1 + ^,^0, 
u!3,o,o,oSo + (^1, 53,0,0,0^0 + ^1 WAflfifl^o + ^1, ^0, 55,0,0,0^0 + ^i}- By giving to the parameters 
Wififij, 5j,o,o,i certain values we can take concrete expressions for the system. For instance, 
if we put wifififi = 16, wi,o,o,i = -10, ^^;2,o,o,o = 1, ^^2,0,0,1 = 1, «^3,o,o,o = -1, «^4,o,o,o = -§, 
5i. 0,0.0 = 2, S3,o,o,o = ~3, S6,o,o.o = ~1 we can eliminate the remainder R^^ and F takes 
the form: F = + 61 - ei-doeQ * [165q - lO^i + 62, 2eo + ei]+ 25oeo * [^0 + ^i + ^2, ^o]- 
lO^oCo * [So + Si, -3eo + ei]+ 52Soeo * [-^Sq + Si, eo]+6So€o * [Sq, -Sq + ei]. 
As a main application of the above methodology, we examine the model-matching problem 
for those systems. This specific problem concerns the discovery of linear feedback connec- 



choice of such a desired system can turn the problem into one that concerns tracking. In 
order to clarify our way of facing up the problem, let us recall some well known ideas from 
the theory of linear systems. Let us have the linear input - output system: 

«oy(^) + ciiy{t - 1) H h aky{t - k) = biu{t - 1) + b2u{t - 2) H h bmu{t - m) 

we want to find a feedback connection so that the closed-loop system behaves like a given 
desired linear system. By means of the operator q, where q~^y{t) = y{t— 1), we rewrite the 
original system and the desired system algebraically as: Ay{t) = Bu(t) , Ady*{t) = B^udt), 
where A, B, A^, B^, proper polynomials of the single variable q~^. We define the feedback- 
law upon request to be of the form Ru — Tuc — Sy, where R, T, S polynomials of q"^ to be 
determined. Using the fact that in the linear case, the usual product among polynomials 
corresponds to the cascade connection of systems (i.e. if = Gu{t) and u{t) = Fv{t) 
then y{t) = GFv{t)) we take the closed-loop system {RA + BS)y{t) = BTudt). Assuming 
that y = y* we obtain the polynomials R, S, T by solving the equations RA + BS = A^, 
BT^Bd, [21]. 

The question, which naturally arises, is how can we extend the above approach in the case 
of nonlinear discrete input-output system. We firstly faced this problem in the paper [11], 
where we introduced the 5-operator in order to replace the q operator and the star-product 
instead of the common product, (therefore the star-product corresponds to the cascade 
connection of nonlinear input-output systems). By means of those tools we managed to 
solve the problem only for input-output systems with nonlinear output, linear input and 
without cross-products between input and output signals. Extension of those results can be 
found at [12], [14], [15], [25]. The aim of the current paper is to present a solution of the model 
matching problem for a general class of nonlinear input-output systems with nonlinearities 
in the output, in the input, and with cross-products. Essentially, what we are doing here 
is to reduce the feedback designing questions to the linear factors Lg, Mg, appeared at the 
original system. Indeed, if a feedback-law u{t) = Sy{t) is available, then the system (1) 
is transformed to a closed-loop form which is nothing else than a difference equation with 
respect to y{t). the core of our approach is to find conditions which guarantee that all 
the terms of the said difference equation have a common linear factor. To achieve that we 
use 5e-operators and some suitable algorithms to express (1) as (2). Then, we throw away 
the remainders R^^,Re^h, R§ and we seek for those values of the free parameters and the 
quantity S, which provide us with linear polynomials Lg, Mg, which are not prime among 
each other. Let us denote their common factor by $. If the desired system coincides with 
$ or it has common factoe with it, then the previous feedback- law solves the problem for 
proper initial conditions. Our empasis is on the fact that the original nonlinear design has 
been reduced to numerous operations among linear polynomials, which can be carried out 
by menas of well known techniques of algebraic geometry. 

Moreover, our approach in the present paper will be purely formal and we will not address 
any issues that pertain to internal stability. We take the main advantages of our approach 
to be the following: (a) that it is a symbohc algorithmic orientation that allows a direct 
computer implementation, (b) that it allows the emergence of a set of solutions for each 
of the model matching problem above, as opposed to a single solution. This is of vital 
importance, since it allows us to choose between available solutions in a way that is informed 
by further constraints. There is strong evidence that this approach, if successful, may 



or optimal controllers design theory. A consideration of these problems, along with an 
assessment of the success of current methods on the basis of specific examples, will be the 
subject matter of a future paper. 

The present paper is divided into two parts. In Part I we present the algebraic frame- 
work and the algorithms, in Part II we are dealing with the model-matching problem. 
Throughout the paper, N, Z, Z+, R will denote the sets of natural, integers, positive inte- 
gers, rational and real numbers, respectively. 

2 The ^-operators and ^-polynomials 

In this section we recall the notions of 5-operators and ^-polynomials. They have been 
appeared and studied in [11], [12], [18]. The scope of their presentation here is to constitute 
an introductory framework for the consideration of the 5e-operators and 5e-polynomials, 
which are going to be adopted at the next section. Let y{t) be a real causal sequence, 
that is a function y : Z ^ R with the property y{t) = 0, ior t < 0. We denote the set of 
such sequences by -F. It is originating from the sampling of continuous functions. Let i 
be an ordered element of the set Z"*"", in other words, i = (zi, ^2, . . . , in) is a vector with n 
positive integers as components, placed in an ascending way («i < ^2 < ■ • ■ < in)- We call 
i a multi-index with dimension n. The ^-operator : F — > F is defined as follows: 

SiVit) = Si^Si^ ■ ■ ■ 5i„y{t) = y{t - ii)y{t -12)- •• y{t - in) 

If in the multi-index i certain components are equal among each other, then we can use 
an alternate notation. Indeed, let ii = ^2 = • • • = = j, ia+i — ia+2 — • • ■ — ia+b — k, 
. . .,in-c = in-c+1 = ' ' ' — in = ^1 then the operator 5; can be rewritten as 5-i = 5"j5\ ■ ■ ■ 51- 
The quantity a + h+ ■ ■ ■ + cis called degree of 5i and it is denoted by deg((5i). If i consists 
from one clement i = (2), then we have the simple ?-shift operator, i.e. diy{t) = y{t — i). 
If at least one of the components of a multi-index is equal to zero then the operator is 
called zero operator. A special case of zero operator is the identical operator Sq, where 
^ouit) = yif)- By convention we define the null operator, 8^ as 5f,y{t) — 1. In other words 
the symbol e corresponds to the " empty " multi- index (). The set of all (5-operators is 
denoted by A. 

Example 2.1 //i = (3), then5,y{t) = 5^y{t) = 5:,y{t) = yit-?.). //i = (1, 1, 1, 2, 3, 3, 4, 4), 
then Siy{t) = (5(i,i,i,2,3,3,4,4)2/(t) = Sl62SlSly{t) = y%t - l)y{t - 2)y\t - ?>)y\t - 4). The 
operators 6oy{t) = y{t), S^S^yit) = ^(0,0,2,2,2) 2/ (^)= y^{t)y^{t — 2) are both zero operators. 

Let I be set of multi-indexes. It may be ordered in a lexicographical way as follows: We say 
that the multi-index i = {ii, 22, • • • , in) is " less " than the multi-index j = (ji, J2, • • • ,jm), 
and we write i -< j if either n < m 01 n = m, and the right-most nonzero entry of 
the vector j - i is positive. For instance, (10) ^ (9,9) ^ (1,1,3,3) -< (1,2,3,3). The 
order among multi-indexes implies an order among ^-operators in a natural way. Indeed, 
we say that 5i -< Si/ if i < i'. Therefore, 5io -< Sl -< Sldl -< 5i52Sl. We equip the 
set of the ^-operators A, with two internal operations: the dot-product and the star- 
product. The dot-product corresponds to the usual product among sequences, while the 
star-product corresponds to the substitution of one sequence by another. Specifically, let 



z{t) = 5,y{t) = y{t-i,)y{t-i2)---y{t-in), w{t) = 5;y{t) = y{t - j^)y{t - j^) ■ ■ -yit - 
im)- We define as dot-product of the operators 5i and 5^ a new operator, denoted by 
5i ■ 5 J, with the property: 5^ ■ Sjy{t) ^z{t)w{t)^ y{t - ii)y{t - i^) ■ ■ ■ y{t - in)y{t - ji)y{t - 
32) • ■ ■ y{i — jm.)- Their star-product is a new operator, denoted by Si*5j, with the property: 
Si*6jy{t) =diw{t)= w{t-ii)w{t-i2) ■ ■ ■ w{t-in) =y{t-ji-ii)y{t-j2-k) ■ ■ ■ y{t-jm-h) ■ ■ ■ 
y{t — ji — in)y{t — j2 ~ ^n) ■ ■ ■ y{^- ^ jrn ~ ^n) ■ In Order to give compact formulas for the 
star and dot-products, we need the following operators among their multi-indexes. Given 
two multi-indexes i = (^1, ^2, • • • , ijt) and j = (ji, J2, • • • , Ja), the new multi-index i©j 
is defined just juxtaposing j after i. Explicitly: i©j = (^1, ii, ^2, ^3, ^2, • • • , 4, Ja) where 
^1 ^ ii ^ ^2 < "^3 < ^2 < ■ ■ ■ < ^fe < 3\- We define the pointwise sum j+z as follows: Let 
j = (ji , ^2 , • • • , jm) be a multi-index and i an integer, then j-j-i = (ji -|- i, J2 + ^, • • • , jm + i) ■ 
Using those notations we get the following formulaes: 

Proposition 2.1 The following properties are valid: (1) 6i ■ 5^ — (2) 5i* — • 

Further properties, concerning the dot and the star products, can be found in [14], [18]. 
We can also define an external operation, named addition, as follows: (^i + dj)y{t) = 
Siy{t) + It can be easily proved that * (^i + Sj) ^ (5k * Si) + (5k * 5j) and thus the 

set (A, +, *) is not a ring. 

Expressions of the form A — J2n=o SiGin '^i'^i called ^-polynomials, where by In we 
denote the set of multi-indexes with dimension n. By convention Iq = {Sg}- For each 
polynomial A we define d{A) as follows: d{A) = min{min(2i, Z2, . . . , in), i = (^i, ^2, ■ ■ ■ , in) ^ 
Z^" such that ai 7^ 0, for n = 1,2, . . . , k}. We define as degree of A and we denote it by 
deg{A) the maximum deg{d) appeared in A. The maximum term of a non-linear polynomial 
A, denoted by max{A), is its largest term, accordingly to the lexicographical order defined 
in the previous section. An expression of the form J2i£z (^i^i is called a linear polynomial. 
Two 5-polynomials A = Yln=o^iein'^i^i ^^'^ ^ — Sm=o SjeJm ^j'^j equal ii 1/ = /i, 
In = Jm, n = 0,1, ... ,u and ai = by Their dot-product is defined as follows: A ■ B — 
I]n=o 2I)m=o SiGin SjGJm '^i^j^i ' '^j- '^^^^ dot-product is nothing else than the classical 
product among polynomials with many variables. In order to define their star-product, 
we proceed as follows: Let us consider the 5-polynomials A and B as functions, trans- 
forming sequences to sequences, accordingly to the rules: A : F ^ F,w{t) Aw{t), 
B : F F, y{t) —>■ By{t), then, the star-product A* B is defined as: 

A*B:F^F , y{t) ^ A* By{t) ^ Ao By{t) ^ A{By{t)) 

There are certain formulaes for the calculation of the star-product. A first one appeared 
in [11], another one is provided in [18]. For the sake of completeness of our presentation 
we shall present the first one, which is the more frequently used. 

Proposition 2.2 Let us suppose that we have the 8 -polynomials A = YZi=o ^i={h,i2,...,i„)&i„ '^i^i? 
B — S(Ji=oSjeJm ^j'^j^ their star-product is given by the relation: 

V 

"=0 i={il,i2,..;in)e1n (jl j2,.--Jn)e(UmJm)" 



Proposition 2.3 (1) * A = En=o EiGi„ ^ ^ Z (2) 6i * A = {6i, * A) ■ {5i, * 

A)--- *A),i = . . .,in) (3) 5,*{A-B) = (di * A) ■ {5, * B) 

The following propositions, can be found in [14]. 

Proposition 2.4 {l)[A + B] * C = A * C + B * C , {2)C *[A + B] = C*A + C*Biff 
C linear, (3)A * B ^ B * A, {A){A * B) * T = A * {B * T), {5)d{A * B) = d{A) + d{B), 
{6)deg{A * B) = deg{A) ■ deg{B) {7){A ■ B) * C = {A * C) ■ {B * C). 

The following property is very useful. It claims that working with linear ^-polynomials and 
the star-product it is like working with polynomials of a single variable and the classical 
product among them. 

Proposition 2.5 Leth be the set of linear 5 -polynomials, then the set (L,*,-|-) is a com- 
mutative ring and it is isomorphic to the ring {R[x\, - ,-{-), where R[x\ is the set of real 
polynomials o a single variable and ■ the operation of the polynomial product. 

Proof: That (L, *, -|-) is a ring comes as a straightforward result of the proposition 2.4. 
For any linear 5-polynomial M — Yli=o frii^i we define the map 

k 

^:(L,*,+)^(i?[x],-,+) , if{M)^Y.^i^' 

1=0 

It can be easily proved that if{M-\-N) = if{M) +ip{N), N another linear 5-polynomail and 
that (fi is one-to-one and onto. Let now A'" = J2'j=o ^j^j-: then M*N — J2i=o J2j=o miUjdi+j. 
This means that ip{M*N) = ELo EJ=o minjx'+^ = Ee±o (ELo msne-s)x^ = (Eto mx') ■ 
{J2'j=oiT'jX^) — <f{M) • <f{N). This relation ensures that is an isomorphism and the 
theorem has been proved. 

3 The (^e-operators and the 5e-polynomials 

The ^-operators defined above, act on single sequences. This concept can help us to 
describe the delays appeared either in the input or in the output signal only. Nevertheless, 
a problem arises when we have to deal with cross-products. Specifically, let us suppose 
for instance, that we have the product y'^{t — l)u^{t — 2), then using 5-operators we get 
5ly{t)62u{t). Hereafter, it is not clear how we shall handle this expression. The main 
difficulty is that we do not know which operator acts where. To pass over this obstacle we 
have to modify the notions of (5-operators and (5-polynomials in a proper way. Actually, we 
introduce the (5e-operator which acts instead of to a single sequence to a pair of sequences. 
This initially appeared in [10] and has been examined in [18]. Let y{t),u{t) be two causal, 
real sequences, defined over the set of positive integers. Let i,j be integers. We define the 
Si X 5j-operator as an operator acting on the pair [y{t),u{t)] (actually, we had to denote 
this pair by {y{t),u{t)), but we use brackets in order to avoid so many parenthesis), as 
follows: 

SiXSj-.FxF^F-.SiX 5j[y{t),u{t)] = y{t - i)u{t - j) 

Let i = (ii, ^2, . . . , im), j = (ji) i2, • • • , in) be multi-indexes. The operator 5i x 5j : F x F i— > 
F is defined as: 



This means that the operator 6i acts exclusively on " outputs " and Sj exclusively on " 
inputs ". Sometimes, for the sake of appearance, the following notation may be more 
convenient: 5i x Sj — SiSy Therefore, the e-operator is just an operator with poperties 
identical similar to the properties of 5-operator, except that it acts only on the second 
sequence (input). Wc call these operators (5e-operators and we denote their set by V. 
Obviously, 6i x 6e[y(t),u{t)] = Siy{t), 6e x 6j{y,u) = ej?^(t). A 5e-opcrator, Sie^ with the 
property that ii = ji = 0, (in other words, the lowest delays of the 6 and e-parts of ^lej 
are zero), is called a zero 5e-operator. A special case of a zero 5e-operator is the operator 
5oeo, with 5oeo[y{t),u{t)] = y{t)u{t). 

Example 3.1 Let i = (0,1,1,2), j = (1,1,3,3,3), r = (0,0,1,2) and s = (0,1,1,2), 
then 5iej[y{t),u{t)] = 5(o,i,i,2) ^{i,i,ifi,2.)[y{t),u{t)] = 5o5i5i52 exexeie^ei[y{t),u{t)] = 5q5\52 
e?e3[y(t),M(t)] =y(t),/(t_l)^(t_2)«2^t-l)M3(t_3)^;^rf/,ermore5ree[?/(t),n(t)] = 5,y{t) = 
SlSiS2y{t) = y^{t)y{t - l)y{t - 2) and Sees[y{t),u{t)] = esu{t) = eoeje2u{t) = u{t)u^{t - 
l)u{t-2) 

The order among multi-indexes implies an order among 5e-operators in a natural way. 
Indeed, we say that ^iej ^ 5iiey if either i < i' or i = i' and j < j'. We equip the set 
of the (5e-operators with two internal operations: the dot-product and the star-product. 
Let 5i X 5j = 5iej, 5ii x 5y = di'ey be two 5e-operators. As their dot-product, denoted by 
^iej • 5i'ey, we define the operator: 

{5i ■ 5i') X {6j ■ 6y) = {Si ■ 6i'){ej ■ ey) = ^ie^ejej' 

It can be easily proved that the dot-product corresponds to the usual product among 
sequences. 

Example 3.2 Let i = (0,1), j = (1,1,2), i' = (0,0,1) and j' = (0,0). Then, Sie^ 
■5i,ey[y{t),u{t)] = 6oS^4f2 dld,el[y{t),u{t)] = 6l6lelele2[y{t),u{t)] = y^t)y{t-l)u\t)u^{t- 
l)u(t — 2). Since i ® i' = (0, 0, 0, 1, 1) and j ® j' = (0, 0, 1, 1, 2) we can easily verify the 
relevant formula. 

As we mentioned earlier, when we were dealing with 5-operators, the star- product corre- 
sponded to the substitution of one sequence into another. In the case of 5e-operators the 
star-product corresponds to the substitution of two sequences with cross-products into the 
"y-part" and the "w-part" respectively, of a given sequence. To clarify this concept, let 
us see an example. Let z{t) = y{t — l)y{t — 2)u^{t — 1)= 5i52e\[y{t),u{t)\ be a sequence, 
involving cross-products of the delays of the sequences y{t) and u{t). Moreover, we suppose 
that each of them is a cross-product of the delays of two other sequences w{t), v{t): y{t) = 
w^{t-2)v{t-3) = 6l€3[w{t),v{t)], u{t) = w{t-l)v{t-l)v{t-2) = 6i€ie2[w{t),v{t)]. By exe- 
cuting the substitution we get: z{t) = w^{t-2)w^{t-3)w'^{t-4:)v{t-2)vlt-3)v{t-4:)v{t-5) 
— ^2<^3<^4e2e3e4e5[w(t), t'(t)], this latter 5e-operator is the star-product of the operator 
5iS2e1 and the pair (62^3, 6iei€2)- Let us present the whole subject formally. Wc have 
the 5e-operators ^iej, Si-^ey, S^^ey. Let [w{t),v{t)] G F x F be an arbitrary given pair 
of causal sequences. We define the maps: $1 : [w{t),v{t)] y[t) = 6i-^ey[w{t),v{t)], 
$2 : Ht),v{t)] ^ u{t) = 6i,ey[w{t),v{t)], $ : [y{t),u{t)] ^ ^(t) = 5iej[y(t), «(t)]. The 
star-product of the operator ^ie; and the pair of operators [(^iieji) ^^12^12] the (5e-operator 



SiiCs, which transfers the pair [w{t),v{t)] to the sequence z{t). In other words, it corre- 
sponds to the map: 

$ o ($1, $2) : [w{t),v{t)] ^ z{t) = Si,e^[w{t),v{t)] 

Alternatively, the star-product is denoted by 5i€j * [^ijej^, (Jijejj]. Of course, by using the 

null operators 6e, wc can include the star-product among simple ^-operators and Se- 
operators into the above definition. For instance: Sie^ * [(^i^ej^, (5eee] = Si * ^iiCjj. Let us 
equip the set of 5e-operators with an external operation, named addition, as follows: 

[S^ej, + 5i^ej^][y{t),u{t)] = Si,ej,[y{t),u{t)] + 5i^ej^][y{t),u{t)] 

The following properties are valid for all the above operations. Their proofs are either 
obvious from the definitions and therefore omitted, or they can be found in [10], [18]. 

Proposition 3.1 (a) The star-product is well defined. 

(b) SiCj * [Si^ej^ySi^ej^] = ^i^+j^^i^^j^ ' ' ' ^h+ir, ^ji+ji^ji+i2 ' ' ' ^h+iu '^i2+ii^i2+i2 ' ' ' 
K+jn %+h%+j2 ■ ■ ■ ^h+ju where i = (ii, is, ... , in) and j = (ji, ^2, • • • , jn) 

(c) 5i * (^i.ejj = (5i * X 5jJ = {5i * (5iJ X {5i * (5jJ = {S-, * 5-J{€i * ejj 

(d) SiSj * [Si.e^, + ^i^ej^j ^ SiSj * Si.Sj^ + S-.e^ * di.e^. 

Example 3.3 LeH = (0, 1), j = (1, 1, 2), ii = (0, 0, 1), ji = (1, 1), is = (0, 1), h = (2), 
then 6iej * [Si^tj^, Si^tj^] = 5o^ie?e2 * [SlSief, 606162]- By setting y{t) = 6l6iel[w{t),v{t)] = 
w^{t)w{t - l)v'^{t - 1), u{t) = 5Q5ie2[w{t),v{t)] = w{t)w{t - l)v{t - 2) and z{t) = 
5Q5i(^\^2[y{t),u{t)] = y{t)y{t-l)u^{t-l)u{t-2), we finally get z{t) = w^{t)w^{t-l)w'^{t- 
2)w{t-3) v\t-l)v\t-2)v^t^3)v{t-A) = 6l6l6^^6s €l4el€i[w{t),v{t)] and thus 5o5ie?e2* 
[6l6 lel, 606 ie2] = '^o'^i^2^3eie2^3^4- Using part (h) of proposition (3.1) we get: (^(o,i)e(i,i.2)* 

[^(0,0,1)^(1,1)' '^(0,1)^2] = (^(0+0,0+0,1+0)^(0+1,0+1,1+1) ^(1+0,1+0)6(1+1, 1+1)5(0+1, 1+1)^0+1, l+l)(^(0+2,l+2) 

c-'2 r4 r 2 2 2 
e2+ie2+ie2+2 = OQOi020ztit2t^ti 

The 5e-polynomials are straightforward extensions of the 5e-opcrators that are necessary 
for the description of discrete non-linear systems with products among input and output 
sequences. These have been introduced and studied in [10], [18]. We will go over the main 
points, adding some new results that will clarify the ideas in the present paper. Let I„, 
be sets of multi-indexes with dimensions n and m respectively. By convention we set 
lo = Jo = {e}. An expression of the form: 

n=0 m=0 (i,j)eln X Jru 

where Cy are real numbers, is called a 5e-polynomial. Whenever a 5e-polynomial acts on a 
pair of sequences, produces a polynomial of delays of those sequences, including also and 
products between delays of the first and second sequence, (the so-called cross-products). 
By means of the null index e, we can decompose a 5e-polynomial into its pure 5-part, e-part 
and Je-part. Indeed, 

/4 r = V V r.....6..f.. = V V r.:.6;f. = V V r:6: 



is the pure 5-part of A and correspondingly 

"^=0 (iJ)eloxJm "i=OjeJr„ m=0j6Jm 

the pure e-part. Actually, it is a 5-polynomial acting only to inputs. The quantity A^^ = 
A — Ag — Ae is the pure 5e-part. We shall call expressions of the form: 

J2 CijSiej = CiSi + J2 cjej (3) 

(i,i)G(iixJo)u(ioxJi) ieii jeJi 

{ij)G(IiXJo) iGli (ij)G(IoXJi) iGJi 

linear de, d and e-polynomials, respectively. The term, which according to the order defined 
previously is ordered highly among the terms of A, is called the maximum term of A. By 
d{A) we denote the minimum delay of yl, in other words d{A) = min(min(i), min(j)), (i, j) G 
In X Jmi n — 0, . . . ^u, m — 0, . . . , fj,. The largest of the numbers deg{i + j), is called the 
degree of A, denoted also by deg{A). The equality of two 5e-polynomials is defined with 
the same manner with the (5-polynomials. 

Example 3.4 For the polynomial: A^ 5SQ + 6di + dl-52Sl + 2ei-3e2 + el-26iei + 36leiel 
we have Iq = {e}, Ii = {0,1}, h = {(1,1)}, I3 = {(2,3,3)}, Jq = {0}, Ji = {1,2}, 

J2 = {(1,1)}; J3 = {(1,2,2)} and Coe = 5, Cie = 6, C(i,i)e = 1, C(2,3,3)e = "1; Cgi = 2, 

Ce2 = —3, Ce(i 1) = 1, cu = —2, C(i^i)(i,2.2) = 3, Cy = for all the other cases. d{A) = 
and deg{A) = deg(5^) + deg(eie2) = 6. When this polynomial acts to a pair of sequences, 
we get: A[y{t),u{t)] = 5y{t) + 6y{t - 1) + y^{t - l)-y{t - 2)y{t - Zf + 2u{t - 1)- Zu{t - 
2) + u'^{t - l)-2y{t - l)u{t - 1) + ?>y\t - l)u{t - l)u^{t - 2). 

Now, working analogously with the case of (5-polynomials we are going to define the two 
internal operations, the dot and the star-product for 5e-polynomials. Let 

^=Y.Y. ai^iA.e^i and S=XX ^ bi2h^i2^h 

n=0 m=0 (ii ji)el„,i X J^,i n=0 m=0 (is J2)eln,2 x J^,2 

be two de-polynomials. Their dot-product, denoted hj A - B, is the de-polynomial, defined 
as: 

max(i/i,^'2) max(/ii,/i2) 

^•^= E E E ci/i^j 

n=0 m=0 (ij)6l„xJm 

with 5i — • 5i2, ej = eji • ejg and cy = Oi^j^ • h^j^. This dot-product corresponds also to 
the usual product among polynomials of many variables. The star-product is also defined 
here as a composition of maps. Indeed, let B, C, A be de-polynomials. We define the maps: 
B : Fx F ^ F,[w{t),v{t)] y{t) = B[w{t),v{t)], C : F x F ^ F,[w{t),v{t)] u{t) = 
C[w{t),v{t)], A: F X F ^ F,[y{t),u{t)] z{t) = A[y{t),u{t)] The (5e-polynomial which 
corresponds to the composition: 



/I „ r D r^i . TP TP .TP 



is called the star-product of the polynomials A, B,C and it is denoted hj A * [B , C]. The 
above definition can be either reduced to some special cases, when B is only a 5-polynomial 
and C an e-one, or extended to more general situations, so that substitutions of a pair of 
(5e-polynomials into a pair of 5e-polynomials, to be described. Details can be found at 
[18]. A question, which naturally arises, is how can we calculate the star-product of Se- 
polynomials. There arc certain formulas, [10], [18]. The most general one is presenting in 
[18], where we can calculate the star-product between pairs of 5e-polynomials, by means 
of tensor-products. Nevertheless, for the problem we study in the current paper we do not 
need this " huge " relation. The formulas, which are sufficient to cover all the cases we 
meet, are the following: 

Proposition 3.2 Let A and B he two Se-polynomials, with A = Yln=o Sm=o S(i,j)GinxJm '^ij'^i^j' 
then: 

(a) Si*A = j:n=0 Em=0 E(ij)Gl„xJ„ Cij^i+.ej+.; « ^ Z+. 

(b) Si*A = {Si^*A) ■ {6^,*A)---{6i^*A), i = (ii, ^2, • • • , 4)- 

(c) 5iej * [A, B] = {6i * A) ■ (ej * B), where by " ■ " we denote the dot-product. 

( d) Let 5iej be a 5-operator, with i = (zi, ^2, • • • , v)? J = (ib ^2, • • • , i\) and L = Y^i&h ^i^i' 
M — Z^jgji iT^j^j linear S and e-polynomials, Ii,Ji<zZ finite sets of multindices, then 

(5i€j * [L, M] = ^ ^s^Ts^ ■ ■ ■ Ts^mh^mh^ ■ • ■ mh^Sg+iei^+i 
self heJ^ 

where s = (si, S2, ■ ■ ■ , s^^), h = {hi, h2, ■ ■ ■ , hx), If, the cartesian products 

If — Ii X Ii X ■ ■ ■ X Ii, — Ji X Ji X ■ ■ ■ X Ji and s+i, h-\-\ the common vector addi- 

^ V ' ^ V ' 

ip— times X— times 

tion. 

Proof: The proofs of (a), (b) and (c) are straightforward. We use these properties in order 
to show (d). Successively we have: 

SiSi * [L, M] = (S, * L) ■ (ej * M) = (5,, * L) • • • (5,^ * L) ■ (ej, * M) • • • (e,, * M) = 





The given formula is nothing else than a reinstall of the above products in a compact form. 

Example 3.5 (1) h^^^l * [25o + 5i,2eo - ei] = \hh * (2(^0 + ^\)\ ■ Vl * (2eo - ei)] = 
(261 + 62) ■{262 + 63} ■(2e2-e3)- (262-63) = 16el6i62-lQe2e36i62+ 4el6 16 2 + 8el6l- 8626361+ 
26l6l + Sei^i^s- Seaea^i^s + 26I6163+ 4.6I6263 - 4.62636263 + 6I6263 

(2) {26162 + 6I6162) *[6i - 2(52, e^ei] = 25ie2 * [^i - 252,egei] +5^eie2 * [61 - 262,6l6i] 
= [2<5i * {61 - 262)] ■ [62 * egei] +[6l * {61 - 262)] ■ [6162 * e^ei] = 2(^2 - 2^3) • {6^63)+ 

{61 - 262)^ ■ 6\62 ■ e|e3 = 2626l63-4636l63+6l6l6l63-46i626l6l63+46l6l6l63 



4 The Formal Linear Factors 



In this section we present the main tool that is used throughout this paper. Actually, it 

is about detecting linear polynomials which can be considered as " factors " of a given 
(5e-polynomial with respect to the star-product and the dot-product. Specifically, let A be 
a concrete 5e-polynomial we write it as follows: 

^ = E • M^ei, * [Le, Me]) + % , + Re,h + R^e (4) 

e 

where L^/, Mq are linear 5 and e-polynomials, correspondingly, (a detailed explanation of 
the above expression will be gived later), including the linear 5 and e-parts of A and 
^, Re,h,Rse remainders. The novelty is that some coefficients of those polynomials 
are not concrete numbers but undetermined parameters, which can take arbitrary values. 
We call these coefficients formal coefficients and the linear polynomials, formal linear fac- 
tors of A. The cornerstone of the current paper is relied on the fact that by giving to these 
parameters proper values, we can take certain expressions of A, appropriate for solving the 
model matching problem. 

We elucidate here that we get the expression (4) and the sets of formal linear factors by 
means of a symbolic algorithm. This algorithm works as follows. First it decomposes A to 
its S. € and 5e-parts. For each such part it finds the corresponding maximum term and then 
it constructs formal linear polynomials Lg, Me and an operator (^i^ej^, so that the operation 
^i^ejg * [Le^Me] eliminates it. Working recursively we eliminate all the other terms but a 
remainder R. This remainder is factorized with respect to the dot-product as i? = ^q^q • R 
and then we repeat the whole procedure, working with R. To establish certain definitions 
to a more rigorous way we need first to present the said algorithm. 

The Formal Linear Factors (FLF) - Algorithm. 

Input: A (5e-polynomial A with no constant terms. 

Output: The quantity h* , the sets i2, O and the polynomials ^, Re,v, = 0, 1, . . . , /i*, 

STEP 0: Set /i = 0, /: = {}, M = {],C = {), O = {}. 

STEP 1: Decompose A as follows: A = A^ + A^ + A^^, where A^, A^, A^^ are the pure S-part, 
e-parts and Se- parts of A, respectively. 

STEP 2: Set u = 1,0; = 0. REPEAT the following substeps UNTIL = or A^ is a 
polynomial of the variable do only. 

SUBSTEP 2a: CALL SUBROUTINE REMAINDER with inputs A^, C,M,C, O, 
h, u,uj and outputs R, C,M,C, O. 

SUBSTEP 2b: Find the maximum integer 9 > such that R = 6q ■ R and R \s a 
(^-polynomial, not containing a constant term. 

SUBSTEP 2c: Set A^ = R,u = u + 1. 

END OF THE REPEAT 

STEP 3: Rename the last value A^ as R^ ^. 

STEP 4: Set n = 0,a; = 1. REPEAT the following substeps UNTIL = or A^ is a 



SUBSTEP 4a: CALL SUBROUTINE REMAINDER with input A^, C,M,C, 0,h, u, to 
and outputs R, £, M, C, O. 

SUBSTEP 4b: Find the maximum integer p > such that R = eg ■ R and R is an 
e-polynomial, not containing a constant term. 

SUBSTEP 4c: Set Ae = R, uj = uj + 1. 

END OF THE REPEAT 

STEP 5: Rename the last value A^ as Re,h- 

STEP 6: IF A^^ = or A^^ is a polynomial of the variables So and eo only THEN set R§^ ^ = Ag^, 
h* = h and goto the Output ELSE proceed. 

STEP 7: Set u = 0,uj = 0, CALL SUBROUTINE REMAINDER with input A^^,C,M, C, 
O, h,u,u and outputs R,C,M,C,0. 

STEP 8: Find the maximum integers K,a > such that R = SqEq • R and ^ is a (5e-polynomial, 

not containing a constant term. 

STEP 9: Set A = R, h = h+1, goto step 1. 

SUBROUTINE REMAINDER. 

Input: A (5e-polynomial A, the sets C,A4,C,C> and the counters h,u,u}. 
Output: The polynomial R. The sets C,M,C,0. 

REM-STEP 0: Set A = 0, i? = A. 

REM-STEP 1: We decompose R as follows: R = Ri + Rni where Ri,Rni are the 
linear and the non-linear parts of R, respectively. 

REM-STEP 2: IF i?; / THEN IF u = 1 THEN set Lo,«,h = Ri, co,u,io,h = 1, 
= 0U {So}, C = CU {1}, C = CU {Lo,u,h} ELSE set Mo,^,h = Ri. co,n,^,'/. = 1, 
O = O U {So}, C = CLI {1}, M = MU {Mo,uj,h} ELSE proceed. 

REM-STEP 3: Set R = Rni REPEAT the following REM-substeps UNTIL Rx = 

or all the terms of Rx become zero terms (i.e. terms of the form SoSi^ ■ ■ ■ Si^eoej^ ■ ■ ■ ej^, 

< i2 < ■ ■ ■ < in, < j2 < ■ ■ ■ jm)- 

REM-Substep 3a: A = A + 1, Rx-i = R 

REM-Substep 3b: Let CA,„,c^,ft5i;,_„_„_^ej^_^_^_^ = cx,u,u.jAA2 ■ ■ ■ ^ir^3i^'n ' ' " ^J- 
be the maximum non-zero term of R\-\, < ii < 12 < ■ ■ ■ < in, < ji < 
j2 < ■ ■ ■ < jra and c\^u,uj,h its coefficient. (Actually, in the first iteration of 
the algorithm, c\^u,oj,h is always a real number. It then becomes a function of 
the unknown parameters Wij^k,h,Si,j,k,h-) 

REM-Substep 3c: IF Si^^^,^ 7^ Sg THEN we form the formal linear 
(5-polynomial 

Lx,u,h = 'Wx^u,h,0^0 + Wx,u,h,lSl + Wx,u,h,2S2 H 1- <^n 

where wx,u,h,o-, 'wx,u,h,i, • • • are unknown parameters, taking values in R and ii 
is the minimum delay of the (5-part of the maximum term ELSE Lx,u,h = 
REM-Substep 3d: IF ej^^^^ ^ ee THEN we form the formal linear 
e-polynomial 



where SA,a),/i,0! SA,a;,fe,i) • • • ai'e unknown parameters that take values in R 
and ji is the mininnum delay of the e-part of the maximum term ELSE 

MA,a;,/i = ee 

REM-Substep 3e: We execute the operation: 

Rx = Rx-l - C\,u,uj,h^0^i2-h ■ ■ ■ ^in-h^0(^j2-ji ■ ■ ■ ^jm-ji * [Lx,u,h, Mx,oj,h] = 
= Rx-l - CA,«,a;,/i'5iA,„,^,hej^ * [Lx,u,h, Mx,u,,h] 

REM-Substep 3f: For Lx,u,h ^ <5e, Mx,^,h / ee we set £ = £ U {Lx,u,h}, 
M = MU {Mx,u,,h}, C = CU {cx,u,u,,h}. = 0U J 

END OF THE REPEAT 

REM-STEP 4: Rename the last value of Rx as R and go to the Output. 

The following definitions will play a crucial role in the context of the model matching 
problem, below. First of all the sets C and Ai, which are formed from all the linear 6 
and e-polynomials that appear in the above algorithm, are called 'the sets of the Formal 
(5-Linear Factors' and 'the Formal e- Linear Factors of A\ respectively. We must note here 
that the use of the word Formal is a little bit excessive, since all the elements of the sets C 
and Ai are not formal linear polynomials. The linear 6 and e-parts of A, usually denoted 
by A^^,A^i, are a first example. The linear polynomial Ri, appeared in REM-STEP 1 
is another. By C* and M* we name the sets £* = £ - {A^^}, M* ^ M - {Ae,i}, C§ 
is the set of the linear factors of the 5-part of A, in other words — {Lx,u,o} C C (it 
includes A^ j), £, A4 are the sets of the 5 or e-linear factors of A which do not participate 

to any cross-product, i.e. C = {Lx,u,h, u ^ 0}, Ai = {Mx,u},h, 7^ 0}, finally C* = 
C — {^4^^}, Ai = Ai — {Ae^i}. The coefficient cx,u,ijj,h is called the companion coefficient 
of the polynomial Lx,u,h or of the polynomial Mx,uj,h and it is denoted by ccoef{Lx,u,h) or 
ccoef{Mx,ui,h)- By convention, ccoef{LQ^u,h) = ccoe/(Mo,tj,/j) = 1 for any u and h. The set 
of all companion coefficients is denoted by C. The polynomials ^, Re,u, i' — 0,1, . . . ,h* , 
R^^ are called the S,e and remainders, correspondingly and the set O the set of the 
formal (5e-operators. Let ^ C £ U be a set consisting from some 5 or e-linear factors of 
A. By C C we denote the set of the companion coefficients of the elements of A. That 
is = {cx,u,uj,h £ C : cx,u,ui,h is the companion coefficient of an element of A}. The next 
theorems describe the behaviour of the FLF- Algorithm. 

Theorem 4.1 The FLF-Algorithm terminates after a finite number of steps. 

Proof: A glance at the algorithm indicates that essentially, it consists from three repeat- 
loops, appeared at the steps 2,4 and the REM-step 3 and a loop defined between the steps 
1 and 9. The termination of those loops ensure the termination of the algorithm. Let A 
be an input of the Remainder-Subroutine. As we pass from the REM-substep 3e we get: 
^new ^previous ~ CA,u,a;,ft^i;v * [Lx,u,h, -^A,a;,ft]) where Cx,u,uj,h is the current maximum 

non zero term. By means of proposition 3.2 (d), we get: 

^'i\,u,h^iKuj,h * [Lx,u,h-, -^A,u),/i] = 5o5i^_i^ ■ ■ ■ 5i„_iieoej2_ji • • • ej^^j^ * [Lx,u,h, Mx,uj,h\ = 



= 5i^5i^ ■ ■ ■ Si^ej^€j2 ■ ■ ■ + other terms 

Therefore, the above substraction will eliminate this 5;. . e,-, . term. The construction of 
the linear polynomial Lx^uM, Mx^u),h guarantees that all the terms of Knew will be ordered 
lower than ^iCj, thus, by repeating the same procedure for those terms, we shall gradually 
ehminate all the non zero terms of A and the REMAINDER SUBROUTINE will terminate. 
Let R be the output of this subroutine. It will be consist of terms of the form: (^of^keoCh- 
Factorizing it by powers of (5oeo (steps 2b,4b,8) we decrease the degree of R at least by two. 
This nesting ensures the termination of all the loops and consequently the termination of 
the algorithm. 

Theorem 4.2 Let A he a given Se-polynomial and A^ ^, Ai^^ its linear 5 and e-parts. Let 
C,M.,C, O , h* and R^^ R^ ^, ReM, h = 0, . . . , h* be the outputs of the FLF- Algorithm, 
after its application to a given de-polynomial A. Then, 

h* K fu.^-.h 

h=Guj=Ou=0 A=0 

(5) 

where ujl,ul,Vu,uj,h, Ku,h, o'^j^h are positive integers, depending from the values of h,u and 
^, cx,u,u>,h e C, e O, Lx,u,h e ^, Mx,oj,h e M and the following " border 

" conditions hold: ki^h = 0, ko,o = 0, ai^h = 0, (To,o = 0, co,o,uj,h = 1; co,u,o,h = 1; 
co,o,o,h = 0, cx,u,c,,h ^Oforu + uj>l, ^i^^ = Sq, 5i^ ,^ ^ ^ = 5e, ej^ ^ = eo, ej^ ^ = ee, 
for u = l,...,ul, u = a = 0, . . . , Vu,u,,h, h = 0, . . . ,h* - 1 and Lo,i,o = A^^, 

Mo.1,0 = Ae^i- Furthermore, this expression is feasible for any 6e-polynomial A and unique, 
under the assumption that the parameters wx,u,h,k, s\,ui,h,k ^'^^ considered as constants. 

Proof: Using reverse recursive relations we can write the evolution of the FLF- Algorithm 
in a condensed form. Let us denote by the outcome of the loop, defined by steps 1 
and 9, at any specific iteration h, Au,h the outcome of the repeat procedure in step 2, and 
Ei_j^h the outcome at the repeat procedure in step 4. It can be easily proved that the next 
equalities hold: 



Hn = Ai,,+£;i,,+^ c,,o,o,/.5i.,o.o..ej,,o,o,,*[i^A,o,^> Mx,o,h]+So'"'^'el-''^'-H^+, , = 0, 1, . . . , h*-l 

A=l 

gu,h ^ 



A=l 



Eu),h — Mo^i^^h + X] ^>^fi,0J,h^i\,<^,h * Mx,u),h + ^o"'^ ■ Eu,+l,h , ^ = 1) • • • ) 
A=l 

with terminal conditions Hh* = RSe,h*^ = R^j^, E^*+i^h = Re,h, h = 0,1,..., h*, 

^h, gu,h, luj,h are positive integers, depending from the quantities h, u, uj. Obviously A — Hq. 
Expanding this relation in details we take: 

51,0 52,0 

A = ^0,1,0 + X] CA,l,0,0^iA,i,a * Lx,lfi + ■ (Lo,2,0 + X CA,2,0,0^iA,2,o * -^A,2,0 + (^o''° " ( 



' ■ iLo,u*,o + XI CA,«*,o,o<^iA,«s,o * Lx,u*,o + Rs^q) ■■■)+ Mo,i,o + X] ca,o,i,o 

A=l A=l 

72,0 ^ 

^JA.I.O * ^A,1,0 + 4''" ■ (^0,2,0 + X CA,0,2,oej;,,2,o * ^A,2,0 + ^q''" ■ (■ ■ ■ 60"°"''° ■ {Mo,uj*,o + 

X=l 

+ C\,0,'^*ofi^h CO* * ^A,a;*,0 + Re,o)) " " O + X CA,0,0,0'^iA,o,o,o ejA.o.o.o * [-^A,0,0, ^A,0,o] + 
A=l ° A=l 

91,1 92,1 

+ C''eo°'' • (^0,1,1 + X CA,lAl^iA,i,i * ^A,l,l + (^o'"' ■ (-^0,2,1 + J2 CX,2fi,l^h,2.i * ^A,2,l + V'' " ( 
A=l A=l 

K*-i 1 

• • • V" ■'•(^0,^,1+ E CA,«t,0,l'^i,,„» i*^A,ut,l + ^(5,i)) • • O + ^O.l.l + E CAAl,iejA,i,i*^A,l,l + 
A=l A=l 

+eo'''-(Mo,2,i + X CA,o,2,iejA,2,i*^A,2,i+eo'''(- • •eo"'~'''(^o,c^^i+ £ CA,o,<,iejA,<,i*^A,<,i+ 

A=l A=l 
Ul 9l,h*-l 

+Re,l)) ■ • O + X CA,OAl(^i;,,o,o,iejx,o,o,i*[-^A,0,l' ^A,o,i]+C''C''( -^0,W-1+ X Ca,1,0,?i*-1 
A=l A=l 

92, h* -1 



A=l 



<^iA,i,/,*-i * -^A,l,/i*-l + ^0'''' '(-^0,2,/i*-l + X CA,2,0,/»*-l'^iA,2,/,*-i * ^A,2,/i*-l + V'** '( 

9ii*, ,h*-l 

h* —1 



<^0 {Lo,ul,_^,h*-l+ X '^A,w;,_^Aft*-l'^iA,.*, ,h*-i*-^A,w;,_^,/i*-l+i?^^^,_l)) 

A=l 



7i,/i*-i 72,?i*-l 

)+Mo,i,h*-i+ X CAAi,?i*-iejA,i,/.*-i*^A,i,/i*-i+eo''''*"'(Mo,2,ft*-i+ X ^mh*-i^h,2,h*-i* 

A=l A=l 

,* 1 A.* 1 n 1 



A=l 

+ X CA,OA/»*-l^iA,o.o,h*-i^jA,o,o./,*-i * [^A,0,/i*-l, MA,o,ft*-l] + RSe,h*))) •••))•••) (6) 



A=l 



Rewriting (6) in a proper compact manner and renaming the quantities gu,h, Ju,h as Vu,ui,h, 
we take the relation (5) upon proving. Now, let A be an arbitrary 5e-polynomial. If it is 



Lo,i,o = p ^0,1,0 = ^e,/, c\,u,uj,h = h, X 0, u,uj 1 and putting the remainders 
equal to zero. If A is nonlinear but contains only (5oeo-terms, then we put R§^f^* = A 
and c\^u,uj,h — 0. In all the other cases we have a maximum term and therefore the FLF- 
Algorithm can be applied. All the above supports the feasibility of the method. Let 
us now deal with the uniqueness question. We explain first that by the expression " 
are considered as constants " wc mean that wc treat the parameters like to be specific 
number. Now, if A has only zero terms or it is linear then the proof of the theorem is 
trivial. Let us suppose that A contains non-zero terms and that, for the sake of proof, 

V'a^ia^jajV'a e R- with \a = (^l.a, «2,a, • • • , «n,a), ja = {jl,a, j2,a, ■ ■ ■ , jm,a) , IS itS maximum 

term. Let us further suppose, for the sake of simplicity, that the second higher ordered 
term of A is ippSi^ej^,ipp G R with = (ii,^, i2,/3, • • • , in,i3), j/3 = i2,/3, • • • , jm,/?) and 
jm,p = jm,a " 1- As WC pass from the REM-Step 3e for the first time we get: Ri = 

A — Cl,0,0,0^ii,o,o,o^ji,o,o,o*[-^1,0,0; ^^l,0,o] and thus A — Ri + Ci,o,0,0^ii,o,o.o^ji,o,o,o*[-^l,0,0, -^i,o,o]- 

«'i,o,o,i'^i + ■ ■ ■ + Sii,a^ ^1,0,0 = siflfi^eo + 51,0,0,1^1 + • • • + ej,_„. By means of proposi- 
tion (3.2), (d), we can see that Ci,o,o,o^ii,o.o,oeji.o,o,o*[^iAO, ^1,0,0] = Ci,o,o,o'^n,„'^i2,a ' ' ' ^^n.a 

ejl.a^ia.a ■ ■ ■ eW +Clfi,0,0{Sl,0,0,jr^,a-l)K,ak2,a - ■ ■ Si„^a ^h,a^j2,a ' ' ' ^3n.,a-l + Othcr tCrmS = 

ci,o,o,o^iaej„ +ci,o,o,o(si,o,ojv„,a-i)'^i;3ej^-F other terms. By equating the coefficients we take 
that ci,o,o,o = i^a-i hence ci,o,o,o is uniquely determined. Now, the maximum term ap- 
peared into Ri is ^i^ej^- A second visit at the REM-Step 3e will give i?i = i?2 + 
C2,o,o,o^i2,o.(),o^j2.o,o,()* [-^2,0,0) ^2,o,o]- Working exactly as before we shall finally get that C2,o,o,o = 
V'/? — ci,o,o,o('Si,o,o,i,„,a-i)- Since this expression is a function of ci,o,o,o and of some parameters 
Wijhki Sijhk considered as constants, we conclude that C2,o,o,o is also defined uniquely. By in- 
duction, we can then see that all the coefficients are uniquely determined. The polynomials 
R^ ^, Re,h, R^^ fj^* consist only of zero terms. These terms arise either from the polynomial 
A or from the products c\^u,Lo,h^i\ „ u> h^h u^h- The unique determination of the coefficients 
c\^u,u>,h entails the uniqueness of the remainders and, thus, the theorem has been proved. □ 

Remark 4.1 The equation (5) can be written shortly as: 

E '^o^C • i^eK^ie * [Le. M,]) + , + R^,, + R^^ 
e 

with 9 is an abbreviation for the vector: {X,u,uj, h). 

Let A be a 5e-polynomial, >C, Al, its sets of the Formal Linear Factors, and C the set of 
the companion coefficients. The elements of these sets, or of any subset of them, can be 
transformed to specific S or e-linear polynomials or to specific companion coefficients, by 
giving to the parameters wx^u,h,k, s\^u,h,k concrete values, according to a set of substitution 
rules. In this case we say that C,M. and C are evaluated over a set of rules U, thus 
writing: ^ ^ tt ^^"^ ^ rr" ^ore rigorously, let W = be the set of the 



u 



w-parameters, written as a vector and r = {a\^i^h,k) a vector of real numbers, which is in 
one-to-one correspondence with the vector W. We say that these parameters follow the 
rule r, thus writing W — > r, if the following substitutions are valid wx^i^h,k = (^\,i,h,k- The 



set C 



is defined as follows: 



aA-l 



A = 1, . . . , A* ^ and ax positive integers , with W = {wx,u,h,k) r} 



N 



{ C 



ri 



r2 



} 



Let A/^ be a set of rules N — {ri, r2, . . . , r^} then L 

Simultaneously, if S = {sxj^h,k) is the vector of the s-parameters and q = {bxj^h,k) a vector 
of real numbers, which is in one-to-one correspondence with the vector S. We say that 
these parameters follow the rule q, thus writing S — > q, if the following substitutions are 
valid sx,j,h,k = h,j,h,k- We define 



M 



{Mx 



bx-l 

Sx,oj,h,k^k 

k=0 



h = 0, . . . , h* ,u! = 1, . 



X— 1, . . . , A* ;j and bx positive integers , with, S 
Let Q be the set of rules Q = {qi, q2, . . . , q^}, then Ai 



Q 



{Sx,ui,h,k) 
{ M 



q} 

M 



qi 



q2 



M 



}. li w and s-parameters appear together, then we say that they follow the rule 

(r,q), if wx,i,h,k = ax,i,h,k and sxjjiM = bx,j,h,k- Let (r, q) 
{cx,u,cvA'^x,i,h,k,Sx,j,h,k) with (W,S) (r,q)}. Let 



(^r, qj, thus writing (W, S) - 
be a concrete rule, then C 



(r,q) 



\NxQ 



(ra,qf,) 



with 



A^, Q be two sets of rules and N x Q their cartesian product, then C\ 
(ra,qfe) E N X Q}. It is obvious that we can apply the above evaluation procedure not 
only to the entire sets JC, M. and C but to any subset of them, even to just one element. 
For instance, I/A,u,/i|ri denotes the evaluation of the 5-polynomial Lx,u,h over the specific 
substitution rule ri. 

Example 4.1 Let us have the 5e-polynomial: 

A^So-2Si + 2S1S2 + A5l + ]-526i - 2828^ + -8l + 3e2 + IGe? - ISeiea - 2ele2+ 



+5e^ + 2eiel 



1 1 11 



1 

16 



451^262 + ^262 + 25i82.e2 — -^^1^2 



(7) 



In order to clarify the notion of Linear Factors, we shall follow the FLF-Algorithm step by 
step. We shall begin by working with the 5-part of A, A^ = 5q — 2b\ + 28\82 + 4(52+|'^253 — 
2828 ^^\8\. Thelm ear part of Ax is 8q — 28i and hence i^o,i,o = ~ 28i. This polynomial, 
which is not formal, is the first member of the set C. Now, the maximum term of the nonlin- 
ear part of Ag is ^8l and the REM-Substep 3e will give: R — A^ — \8^ee * [-^1,1,0; -^1,0,0] 
with Li^i^o = «'i,i,o,o'^o+ ^i^i,i,o,i'^i+ «'i,i,o,2<^2+ w'l, 1,0,3^3 + ^4; M,o,o = ^e, tMs Li_i,o is 
the next member of C and ccoef{Li^ifi) = —j, O = O U {S^ee}- Now, the higher or- 
dered term of the obtained polynomial R, is — ^^'^^"'^ 8384 and the REM-Substep 3e will give 

Rnew = Rold + "^^'2""' ^O^l(^e * [-^2,1,0, ^2,0,o] ^^^^^ -^2,1,0 = W^2,l,0,0'^0+ W^2,1A1'^1+ W^2,l,0,2<^2+ 83 

and M2fifi — e^- Working similarly we shall finally get the neoct epression for A^: 

Ax ^80- 28i + ^8l * (iWi^i^co^o + «^1,1,0,1<^1 + «^l,l,0,2<^2 + «^l,l,0,3<^3 + 84) - ^^'^'°'^ (5o(5i* 



/ _ Uii 1 n 9 U)i 1 n -iWo 1 n 9 \ . 



*{w3^1fl^o6o+Ws^lfi^lSl+62)+C45o64*{w4,l,0,oSo+Si)+C56l*{w5,l,0,oSo+W5,lfi^^^ 
+C65o(^l*(w6,l,0,0^0+W6,lAl^l + <^2) + C7(Jo(^2*(w7,i,o,o5o + 5l) + C8^0*(^8,lAO<^0+W8,l,^^ 
+Cg5Q5i * (^9,1,0,0(^0 + (^l) + CiqSI * (1010,1,0,0^0 + 5i) + R 

where R contains only zero terms, particularly R — To^q + ri^of^i + r25o<^2 + T^iS^S^. The 

coefficients C4, C5, Ce, Cy, cg, cg, Cio,ro, ri, r2, are polynomial functions of the wx^jis pa- 
rameters. They are not presented explicitly due to their large size (the coefficient Cio, 
for instance, contains 77 monomials) . We factorize R by do and we get: R = 6(,R — 
^o(''o^o + ''i<^i + ''2<^2 + ''3<^3) o-nd we repeat the procedure with — ro5o + ''i<^i + ''2<^2 + ''3^3- 
Since it is linear the SUBROUTINE REMAINDER will terminate at REM-STEP 2, giv- 
ing 1/0,2,0 = '''o'^o + ^^i^i + ''2'^2 + ^3^3 and ccoc/ (I/o,2,o) = I, O = O U {So}. All the above 
linear polynomials are elements of C We are now in Step 4, working with the e-part of A, 
Ae = 3e2 + 16ei — 18eie2 — 2e\e2+ 5e| + 2eie\—\e\ + 6163—616263 + |e|63. Initially, we get 
M),i,o = 362 and for the nonlinear part of A^ we shall finally get: 

Ae,nl — ^^0^1 * (*1, 1,0,0^0 + Sl, 1,0,161 + 62) + ^ — 1 ^2^) * (^2,1,0,0^0 + 61) + 

, /1 1,0,1 , , -^l, 1,0,1-52,1,0,0 \ 2 / , \ , f ^ -51,1,0,1 \ 3 
+ (^1 + S2,l,0,0 + " j eo€2 * (S3,l,0,060 + 61) + [-- j 6^* 

*(s4,i,o,oeo + 54,1,0,161 + 62) + c^eoej * (55,1,0,060 + 61) + CQelei * (s6,i,o,o6o + 61) + 
+C7el * (57,1,0,060 + 61) + 56o * (s8,i,o,oeo + S8,i,o,iei + 62) + (-18 - 10s8,i,o,i)6o6i* 

*('59,l,0,0C0 + Cl) + C1060 * (Sio,l,0,oCo + 61) + i? 

The coefficients C5,C6,C7,Cio have large size and therefore are not written in details. The 

quantity R is R = r2,36o6263+ ri,26o6i62+ ro,362e3+ ro,i6^6i+ r2,26o6|+ ri,36oei63+ ro,o6|j + 
^1,16061 +'"0,26062 +^"26062 + ^"16061 +^o6o- Factorizing R by 60 we get R = eo- R and then we 
repeat the substeps of step 4 once more, taking uj = 2, Mo,2,o = ''"060 + ^"161 + ^262 and R = 
^2,36263+ ri,26i62+ ro,36o63+ ro,i6o6i+ r2,26i+ ri,36i63+ ro,o6o+ri,i6i +ro,26o62. The running 
of the substeps of REM-STEP 3 will finally give R — ro,i6o6i * (si,2,o,o6o + ^1,2,0,161 + 62)+ 

^'0,26062 * (S2,2,0,060 + 61) +ro,o6o * (S3,2,0,060 + S3,2,0,l6l + 62) +ro,i6o6i * (s4,2,0,o6o + 61) +^0,060 * 

(55,2,0,060 + 61) + i? With R = foel + fieo6i + r26o62 + r3eo63 and 7\o,ro,i,ro,2, ro.h.h.h 
are polynomials of all the above parameters. Visiting SUB-Step 4b once more we shall take 
R — Eq- {fo^o + ^1^1 + ^262 + ^363). Now, u) — ?> and since R is linear the SUBROUTINE 
REMAINDER will terminate by giving Mo,3,o = ro6o + ri6i + f2e2 + 7^363. All the above 
appeared linear polynomials are elements of M. . We go now to step 5, working with the 
6e-part of A, A^^ = Si62ei + l6lei+ ^SiSsei -^5361-451(5262 + 52^2+ 2515362-^5362. We 
have u = 0,u = and successive applications of the substeps of REM-STEP 3, will finally 
give the following expression for A^^: 

A§^ = -^^leo * [w^i,o,o,o5o + Wi,o,o,i5i + wi,o,o,252 + 53, si,o,o,oeo + si,o,o,iei + 62]+ 

+ (~Y^ + ^i^^^ 5o6o * [w2,o,0,o5o + «^2,0,0,l5l + ty2,0,0,252 + 53, S2,0,0,060 + 61] + 



H ^^5o(^ieo*[w4AO,o(^o+W4AO,i'^i+'^2,S4AO,oeo+S4AO,iei+e2]+P5^o5ieo*['i^5,OAO^o+'i^5,o,o,i'^i+ 

+^2, 55,0,0,0^0 + ei] + peSoSieo * [■u;6,o,o,o(^o + ■w^6,o,o,i(^i + ^2, eo]+ 

, /r, , ^1,0,0,1 Sl 0,0,2*4,0,0,1 \ r r r X , X i i 1 i 

+ ( 2 H ^ ^ 1 dod2eo * F7,o,o,ocio + Oi, S7,o,o,oeo + S7,o,o,iei + £2] + 

+P8(^0^2eo * [lW8,0,0,0^0 + Si, 58,0,0,0^0 + Cl] +^9^0(^260 * [^"9,0,0,0(^0 + (^1, Co] +Pl0(^o^O* 

*[w^io,o,o,o(^o + ^i^io,o,o,i^i + ^2, sio,o,o,oeo + sio,o,o,iei + €2] +Pii5o^o * [^^^11,0,0,0^0 + «'ii,o,o,i^i + ^2, 

Sii,o,o,oeo + ei] + ^12*^0^0 * [^^'12,0,0,0^0 + it'i2,o,o,i^i + ^2, eo] + Pis'^o'^ieo * [^13,0,0,0^0 + ^1, ^2] + 

+Pi4(5o5ieo * [«'i4,o,o,o'^o + ^1, si4,o,o,oeo + ei] H h -R 

where R = go^o^o + gi^^c^ieo + q2SoS2eo + gs^o^seo- Since R = 6oeo ■ R = Soeo{qoSo + qiSi + 
12^2 + Q3S3) the algorithm will terminate with i^o,o,i = Qo^o + ^Zi^i + Q2S2 + Qs^s- Summa- 
rizing, the final output of the algorithm is C — {Sq — 25i, uji^i^co^o + it'i,i,o,i<^i + Wi^i_o,2<52 + 

^i^l,l,0,353 + ^4, W2,l,0,0<^0 + «'2,1,0,1'^1 + ^2,1,0,2^2 + ^3, • • •}; M = {862, Sl,l,0,oeo + S^i^o.lCl + 
^2, 52,1,0,0^0 + ^1) • • •}; C> = {5q, ^4, SqSi, ...,60, e^ei, . . . , Sq€q, . . .}, C = — H!liM^ _2 — 
M)!, 1,0,2 I «'i,i,(),3U'2,l,0,2 1-1 SI, 1,0,1 ~ 11 I si,0,0,l n n 1 nnfl 
2 ' 2 ' "^41 C51 • • • ? 4) "J- 2 — ' • • • ' ^5) • • • ~4) ~i6"l 4 — ' • • • )P5iP5) • • •/ f*^" 

-Re,v = Rs „ = 0, h* = l,Rsefi = -^Je,i = 0- Moreover, C§ = C- {ro6o + ri6i + r2S2 + r3S3} , 
jC-{do-26i}, M* ^ M-{3e2}, C =_>C- {t^i,o,o,o(^o + wi,o,o,]A + tyiAp,2<J2 + 53, • • •}; 
M ^ M - {si,o,o,oeo + si,o,o,iei + £2, • • •}; jC* ^ £ - {Sq - 25i}, M* ^ M - {Sea}. Let 
us consider now the rules ri = {wh,u,u>,x = 1}; ^2 = {w^i.i.ca = 0,1^1,1,0,2 = -4,W3,i,o,i = 0, 
'W'4,1,0,0 = 0,1^1,1,0,1 = 0, all the other w-parameters take arbitray values }, qi = {sh^u,Lu,\ — 
1}; c[2 = {si, 0,0,1 = all the other s-parameters take arbitray values }, then £|ri = 
{^0 - 2^1, 50 + ^1 + ^2 + 53 + So + 5i + S2 + 6s, . . .}, £^|r2 = {5o - 25i, u^i, 1,0,0^0 - 

4^2, ^2,1,0,0^0 + «^2,l,0,l^l + «'2,1,0,252 + 53,.. A^\r^ ^ Sq - 25i + * (Wl, 1,0,0^0 " 452) + 

|5o5i*(w6,i,o,o(^o + ^y6,i,o,o(^i + (^2) + -R|r2; A1|qi = {862, + 61 + £2, + Ci, . . .} and finally 

C2,0,0,o|(r2,q2) = ("l^ + ^i^) |(r2,q2) = 

5 The Model Matching Problem 

In this section we present the main apphcation of our methodology. Before we examine 
it, we would like to begin by a short algebraic description of the non-linear discrete input- 
output systems, via the notion of 5e-polynomials. Suppose that we have a system of the 
form (1). By using 5e-operators we can rewrite this as follows: F[y{t),u{t)] — 0, where F 
is a 5e-polynomial. As usual, many times we decompose F to its pure 5, e and 5e-parts, 
i.e. F = F^ + F^ + F^^ and F^, F^ to their linear and nonlinear parts: F^ = F^ ^ + F^ 

= F^i + F^^ni- is already purely nonlinear). The causality and solvability of the 

system is guaranteed by the equality d{F) = d{F^ J. Often, we shift all the delays properly 
so that d{F) — d{F^ ^ — 0, which means that the lower delayed term, i.e. y{t), appears 
in the linear 5-part of the system and thus we can solve (1) with respect to y{t) in a 
direct way. To each nonlinear system of the form (1) we assign a vector of real numbers 
Yo = {yoyUii ■ ■ ■ iVk-i) which gives the so called initial conditions: y(0) = y^, y{l) = yi, 
. . . , y{k — 1) = yk-i, where k is the maximum delay appeared in the output signal. Since 
the signals, involved in (1), are causal, i.e. y(t) = 0, u(t) = 0, for t < 0, for each given 



satisfies (1), for n > k. Let K = J2i=Q Ti6i be a linear 5-polynomial, we say that the 
initial conditions yo satisfy the relation Kyo = if Yli=o TiVk-i-i = 0. Let us have two 
causal nonlinear systems: F[y{t),u{t)] with yo = (yo, yi, • • • , Z/it-i) and F[ip{t),v{t)] with 
V^o = ii^o, i^i, • • • ; V'A-i) and k > X. We say that the two systems " operate " under identical 
initial conditions, if ipo = yo, ipi = yi, . . ., ipx^i = yx-i and yx = i^iX), yx+i = i'iX + 1), 
. . ., yk-i = ip{k — 1). In other words wc must give as initial conditions to the system which 
starts to product outputs later, the corresponding outputs of the other system. The next 
theorem is crucial for the model mathcing problem. Its proof can be found at [9] . 

Theorem 5.1 Let Gy{t) = and ^y*{t) = 0, G, $ a nonlinear and a linear 6 -polynomials 
be two systems without inputs, a nonlinear and a linear one, (actually they are difference 
equations with respect to y{t) and y * (t) ). If G — G * ^, then y(t) — y*{t) whenever the 
two systems operate under identical initial conditions. If G = G * ^ and $ = $ * then 
y{t) = y*{t) whenever the two systems operate under identical initial conditions yo, which 
satisfy the relation $yo = 0. 

The model matching problem for input - output discrete systems has already attracted a 
lot of attention. Though there is already a clear understanding of the linear problem [17], 
[21], the non-linear case is still being examined and is the subject matter of a large number 
of papers (see for instance [11], [25], [22], to mention but a few). 

In this paper we examine the Model Matching Problem (MMP) for non-linear discrete 
systems of the form (1), or F[y{t),u{t)] = 0, F a 5e-polynomial, by using the entire 
algebraic background that has been previously developed. This approach is a continuation 
of research on non- linear discrete systems in [10], [11], [14], [12], [15]. Our aim is to find a 
feedback relation of the form u{t) — Sy{t), S a linear (5-polynomial, so that if it is fed back 
to the original system we get an output that is identical to that of a given desired linear 
system, described by the equation Ady*(t) = 0, a linear ^-polynomial. If this holds, we 
say that the nonlinear systems matches the linear one. 

To achieve a solution of the problem we work as follows: First by means of the FLF- 
Algorithm we obtain all the linear factors of the nonlinear system and then we rewrite 
F as in (5). By substituting the feedback relation u(t) = Sy(t), we eliminate the input 
signal u{t) and we take an expression containing only the output signal y{t), actually we 
have a difference equation with respect to y{t). The majority of the terms, appeared in 
this difference expression contain either the formal 5-linear factors or the formal e-linear 
factors of F or the " free " linear polynomial S. The cornerstone of our approach is 
the following: if we can achieve to find values of the parameters Wijhk, Sijhk and of the 
free linear polynomial S, so that all the terms of the said difference equation (formal and 
no-formal), have a common factor which is a factor of A^, as well, then, these specific 
values of S will constitute a solution of the model matching problem. To accomplish this 
task we examine several cases. In all of them wc eliminate the unneccesary remainders 
and some of the linear factors and then we " ask " from the linear parts of the original 
systems (which do not have parametrical coefficients and thus cannot be influenced by 
the choice of proper values) to have a common factor. This can be achieved by solving 
certain Diophantine equations. Two things are remarkable here, first we do not get a single 
solution but a class of solutions, parameterized via a linear polynomial Q. Proper choice 
of this polynomial can create a causal feedback-law or a law with as much or as few delays 
we like. Secondly, and most important, all the calculations and procedures are devoted to 



in section 2 the star-product of linear 6 or e-polynomials coincides with the usual product 
among univariate polynomials. Hence, we can use all the available results and tools of the 
classical algebraic geometry or combutational algebra, to carry out the manipulations we 
need. Before we present a symbolic algorithm, named the MM-Algorithm, which solves 
the problem completely, we shal describe the G CD- Values Algorithm. This algorithm will 
be used as a subroutine of the main algorithm. What it essentially does is to discover 
those values of the parameters which permit to a set of formal linear polynomials to have a 
common factor. To achieve that it works in two steps: firstly it discovers the values of the 
parameters which make certain companion coefficients no- zero and thus the corresponding 
hnear factors are present in the expression (6) and then it selects among these values, those 
specific ones which allow the corresponding linear polynomials to have a common factor. 

The GCD- Values-Algorithm 

Input: The set a formal linear polynomials J^. 

Output: The set of rules GCDV{J^). 

STEP 0: We set GCDV{J^) = {} and m = card{J^), that is m is the amount of polynomials 
contained into T. 

STEP 1: FOR n = 1 TO m 

Step la: We construct the power set P, consisting from all the subsets of T, with n 
elements, i.e. P = : .4 C ^ with card{A) = n}. 

STEP 2a: REPEAT for all ^ G P. 

STEP 2a-I: Construct the set C-^ of the companion coefficients of the 
elements of A. 

STEP 2a-II: We form the set of rules NZ = {r : r is defined by the substitu- 
tion {wijhk,Sijhk) = {aijhk,bijhk), aijhk,bijhk G R. such that ccoef{J) 0, 
for J G ^ and ccoe/(*) = for ^ e - A}. 

STEP 2a-II: We form the set of rules CF C NZ such that gcd(Ji) 7^ 
5o,Ji e A\nz and we set GCDV{J^) = GCDV{J^) U CF. 

END OF THE REPEAT 

END OF THE FOR 

Remcirk 5.1 1) The calculation of the values of the parameters which ensure thatgcd{Ji) ^ 
So, can be carried out by means of resultans theory [2]. 

2) Obviously, if gcd( Jj) = 6q, for Ji belonging to a specific set of polynomials A, then 
for any set B with A <Z B, we shall have gcd(Jj) = 60, Ji e B, too. Therefore, we can 
avoid to check any set which includes A. This comment can reduce the whole amount of 
calculations, significantly. 

We are now ready to state the main algorithm of the paper. 



The MM-Algorithm 



Output: The sets Fi,F2,F3. 



STEP 1: We decompose the polynomial F into its linear and non-linear parts: F = 

Fs,i + ^^,1 + ^ni and we rename A = F^ ^ and B = Fe^i. 

STEP 2: By means of the FLF-Algorithm we form the sets JC, M, C*, M*, C§, Z, 

M, €*, M* and R^^ = Re,h, h = 0,...,h*, R^^. 

STEP 3: By means of the GCD-Values-Algorithm, whenever it is needed, we form the 
following sets of rules: 

• LV = GCDV{C^) 

m GV = GCDV{C U A7) U GCDV(C U M) 

• MV = GCDV{C* U M*) U GCDViZ* U M*) 

• DR = {{wijhk) = {a-ijhk), aijhk G such that i?j ^ = 0}. 

• RV = {{wijhk,Sijhk) = {aijhk, hjhk), aijhk,bijhk ^ R. such that Rg f^ = Re,h = 
0, h = 0,...,h*, Rg^ = 0} 

STEP 4: IF LV n DR $ THEN put Fi = {$ : $ a common factor of Lj G 
STEP 5: IF RV nGV ^ THEN put F2 = : $ a common factor of Li G 



and Mi e M 

RVnGV 



and Mi e M 



} U 1$ : $ a common factor of Lj G C 

RVnGV^ 



RVnGV 



RVnGV^' 



STEP 6: IF RV fl MV THEN put F3 = {$ : $ a common factor of 

Li e C 



and Mi e M 

RVnMV 



I U 1$ : $ = a common factor of 



and Mi e M* 

RVnMV 



RVnMV^ 



LiE r 

STEP 7: Goto to the Output 

The next theorem includes the feedback construction upon request. 

Theorem 5.2 Let F[y{t),u{t)] = 0, F a de-polynomial, be a nonlinear discrete input- 
output system with cross-products. Let A be the linear output part of F, —B its linear 
input (i.e. F^i — A,Fe,i — —B) and Ady*{t) —0,Afia linear 5 -polynomial, a linear desired 
system. Let Fi,i = 1,2,3 be the outputs of the MM-Algorithm, after its application to the 
polynomial F. We suppose that at least one of the sets Fi,i = 1,2,3 is novoid. Then: 

a) If Ad belongs to some of the sets Fj,z = 1,2,3, then the feedback law u{t) = Sy{t), 
with 

r Q*Ad, tfAdeFi 
S=\ Q, tfAdeF2 
[ Zo + Q*Ad, ifAdEFs 

where {Rq, Zq) is a solution of the Diophantine equation R * A^ + Z * B = A, and Q 
an arbitrary linear 6 -polynomial, gives a closed-loop system which produces an output y{t) 
equal to y*(t), whenever the nominal nonlinear system and the linear desired one, "operate" 



b) If Ad does not belong to any of the sets Fi,i = 1,2,3 but there are subsets Fj C 
Fj, i = 1,2, 3, not all of them void, such that gcd{Ad, $) = $ 7^ Sojor $ G Fj, i = 1,2, 3, 
then, the feedback law u{t) — Sy{t), with 



S 



' g * !>, $ e Fi 

Q, z/<fGF2 

, Zo + Q*$, z/$ e F3 



where (Ro, Zq) is a solution of the Diophantine equation R * ^ -\- Z * B — A, and Q 
an arbitrary linear 5 -polynomial, gives a closed-loop system which produces an output y{t) 
equal to y*{t), whenever the nominal nonlinear system and the linear desired one, "operate" 
under identical initial conditions yo, which satisfy the relation $yo = 0. 

Proof: Let F[y{t),u{ty\ = be the nominal plant, F a 5e-polynomial. We shall follow 
the MM- Algorithm step by step. By means of the FLF- Algorithm and theorem 4.2 we get 
the following expression for F. 

h* K '"u,^,h 

h=GijJ=Ou=0 A=0 

+% + Re,ft + R^^ (8) 

Using the border values, given at 4.2, and the fact that -^^0,1,0 ~ ^5 1 ~ ^1 -^0,1,0 = Fe,i — 
—B, we can split (8) as follows: 

«(*) fu,0,() t^O f(),u),0 

F = ^+Y1 Yl ^o"'°-CA,«,0,0(^iA,„,o.o*^A,^i,0 + %_0~^+I^ Yl eo"'"-CA,0,c^,oejA,o,a.,o*^A,u;,0 + 
u=l A=l ' ui=l A=0 



^0,0,0 K <J ^* 1 ^"■"■'^ 

+Re,o + E CA,o,o,o^iA.o,o,o^JA,o,o.o * [^A,o,o, M,,o,o] + <5o"°'°eo"°'° E Lo,i,H + E E ^o"'" 

A=l h=l \ u=l A=l 



ui=l A=l 

+ E ^0 eo • CA,o,o,?^<^iA,o.o,/.ejA,o,o.h * [L\fl,h, Mx,o,h] + (9) 

A=l / 

The above expression (9), recovers the " output " part, the " input " part and the " cross- 
products " part, contained into the original polynomial F. By substituting (9) and the 
feedback relation u{t) — Sy{t) into the original nonlinear discrete system, we turn it to the 
following difference equation with respect to y{t): 

Mt) + i: E <5o"'°y(i)-CA,„,o,o^iA.„.o,o*^A,«,oyW+%,o^W-5*5?/(i) + E E ^o'"°*Sy{t) 



""0,0,0 

CA,0,0,oOi;^_o,o,o^jA,o,o,o 

A=l 

+5o""%(*)eo"'''° * Sy{t) Y: Lo,^Mt) + E E C''y{t) ■ CA,„,o,/.^i.,.,o,. * L,,^Mt)+ 

h=l \ u=l A=l 

+Rs^^y{t) + Mo,i,h*Sy{t) + E E ^o"'" * "^^/It) • CA,o,c.,/^ej,,o,.,. * Ma,^,/.* +Re,„ * Sy{t) + 

ui=l A=l 

^'o,o,(. K a \ 

+ + E '^o""' ?/(t)6o"''' * Sy{t) *[LA,o,/,y(t),MA,o,ft*5'?/(i)] + 

A=l / 

Rs,*S[y{t),u{t)]^0 (10) 

Where we used the fact that A, Lx,u,h, R§ ^ ^-re polynomials of 5-operators only, and thus 
they act exclusively on the output, whilst B, Mx^u,h, Re,h ^-re polynomials of e-operators 
only, and thus they act exclusively on the input. By means of proposition 3.2 (c) we can 
easily prove that: 

(i) (5iej * [Ly, M * Sy] = {5i * L) ■ [5-^* M * S)y{t), for any linear polynomials L, M, S (the 
operator has changed to 5j, since we are dealing only with y sequences). 

(ii) * Sy{t) = S^y{t) 

(iii) Re,H * Sy{t) = c^A" * Syit) = E,, C0,^S^'y{t). 

(iv) = E(a,/3) C(a,p)^o4*[yit),Sy{t)] = E(a„fl) C{a„3)K Sl^yit) = y{t)] = E(a,/3) C{a,l3)y''{t)[Sy{t)f 

Therefore, after some manipulations, (10) becomes: 
{A-B*S) + J2i: So"-" ■ cx,ufl,oSi.,.,o,o * Lx,ufl + %,o + E E ca,o,.,o5'^-° • (5j.,o,.,o* 

u=l A=l a;=l A=l 

fO,(),0 

*Ma,,,0 * S) + Ec/So-?* + E CA,0,0,o((^i.,o,o,o * ^A,0,o) " (<Jj,,o,o,o * ^A,0,0 * 5) + 

• E ( ^0,1,/. + E E <^0"''^ ■ CA,«,0,.<^i.,.,o,. * L,,^,h + 
/i=l \ u=l A=l 

a;=l A=l 



+ E CA,OA^<^o"" S^-l^' ■ (5i,,„,„,, * La,o,.) • (5j,,o,o,. * Ma,0,. * 5) + E Cia,P)S"oS^ 
A=l / (a,/3) 



t/(i) = 



The left-hand side of (11) is a polynomial of (5-operators only. We denote it by Ff. So, 
(11) is nothing else than the closed-loop system we take from F[y{t),u{t)] = by using 
the feedback-law u{t) — Sy{t), in other words F[y{t), Sy{t)] — Ffy{t) — 0. We are now 
ready to state the main argument of this proof. Since the desired system Aayit) = has 
no input, it is obvious that the closed-loop system will have the same dynamic behaviour 
as the desired one, only if we can find values for the paramctrical coefficients and the 
unknown quantity S, such that the difference equations (11) and Aciy{t) = have common 
solutions. To achieve this, we first need to use the equation (11). By means of theorem 
5.1, we can establish that if we can find a set of values of the parameters so that $ is a 
linear common factor of all the evaluated terms of (11), then any solution of the difference 
equation ^y{t) = 0, is also a solution of the equation (11). If we ensure that this $ is a 
factor of the polynomial Ad, too, or it coincides with it, then we have the desired result. 
Therefore, we have to consider the following three cases, each of which corresponds to one 
of the three steps 4,5,6. 

Case I. Let us suppose that Fi is non void and that $ G Fi. Step 4 of the MM-Algorithm 
implies that LV fl DV ^ 0. This condition means that we can find at least a rule of 
substitutions for the parameters {w\^u,h,k, s\,u),h,k), named q, such that two targets are 
accomplished simultaneously: (i) The remainder ^ is eliminated. This is due to the fact 
that q e DR. (ii) There is a certain amount of polynomials ivA,«,o with ccoe/ (LA,u,o)|q 7^ 
(thus the polynomials La.u.oIq ^-rc not anihilated from Fj|q and hence we can calculate 
their great common divider ) and gcd(LA,M,o|q £ ^§\ci) = ^ ^ ^o- This is due to the fact 
that q e LV and the construction of the GCDV- Algorithm. Now, we set S = Q * ^, Q 
an arbitrary linear polynomial. Since A e £^|q , (A has not parametrical coefficients) =^ 
A = T * for some T and hence A-B*S^T*^-B*Q*^^{T-B*Q)*^^K*^ 
and so $ is a factor of A — B* S, too. Observing that is a factor of all the terms of (11), 
except the terms of the first line and taking into account all the above posed arguments 
and the structure of Ff, we finally conclude that $ is a factor of all the terms of i^/|q, 
in other words -P/|q = * $. Now, if = $ then, by means of theorem 5.1 we get 
that the nonlinear system matches the desired linear one under any initial conditions yo. 
If gcd(Arf,$) = $ ^ (^0 then, by means of theorem 5.1 we get that the nonlinear system 
matches the desired linear one under any initial conditions yo, which satisfy the relation 
$yo = 0. All the above arc valid for any polynomial $ € Fi and hence the first branch of 
the feedback law of theorem 5.2 has been proved. 

Case II. In this, rather extreme case, we want to find a common factor $ of the terms of 
(11), by leaving the polynomial S totally free. Let us suppose, as before, that F2 is non void 
and that $ G F2. Step 5 of the MM-Algorithm implies that RVOGV ^ 0. This condition 
ensures that there is at least one rule q such that all the remainders are eliminated and 
there is either a subset of the set (£ U 7M)|q or a subset of the set {CU A4)\q consisting 
from polynomials P, with no zero companion coefficients and with the property: either 
gcd(P G (£uAT)|q) = $ ^ ^0 or gcd(P G (SuA<)|q) = $ ^ 5o- Since A G C^A = A*$ 
and B G M^B = Si * $, for some polynomials Ai, Bi, so A - B * S= (Ai - Bi* S) * ^, 
which implies that $ is a factor of A — B * S, too. From all the above we conclude that 
$ is a common factor of all the terms of -F/|q, i.e. F/-|q = * $, independently from the 
values of S. Therefore, if A^ — ^ then by means of 5.1 the nonlinear system and the linear 
desired one, will have the same dynamic behaviour for any initial conditions and any value 



means of 5.1 the nonlinear system and the hnear dcsh^cd one will have the same dynamic 
behaviour for any initial conditions yo, which satisfy the equation $yo = and any value 
of S, thus we set again S — Q, Q an arbitrary linear polynomial. All the above are valid 
for any polynomial $ G F2 and hence the second branch of the feedback law of theorem 
5.2 has been proved. 

Case III. We arc dealing now with a case similar to the second one. The essential difference 
is that the linear parts of F, which do not have paramctrical coefficients are not involved 
to the calculation of the common factors. Let us suppose, as before, that F3 is non void 
and that $ G F3. Step 6 of the MM- Algorithm imphes that RVDMV ^ 0. This condition 
ensures that there is at least one rule q such that all the remainders are eliminated and 
there is either a subset of the set (£*U7Vf )|q or a subset of the set (£ U7Vl*)|q, consisting 
from polynomials P, with no zero companion coefficients and with the property: either 
gcd(P G {L* U Al*)|q) = $ ^ 5o or gcd(P G {T U Al*)|q) = $ 7^ ^o- These relations 
mean that $ is a common factor of all the terms of P/|q, but the term A — B^S. We have 
to choose now S properly so that this latter term is a multiplier of $, too. This means 
that A — P * S* = i? * $, for some polynomial R and thus P*S' + i?*<l> = A. This is 
the well-known Diophantine equation among linear polynomials, which has a solution if 
gcd(S, ^)\A. If gcd(i?, $) = 5o the above assumption is true. If gcd(i?, <I>) = l> ^ 5o, then 
we conclude that gcd(yl, P,$) ^ but this implies that gcd(P G (£ u!M)|q) = $ 7^ 5o 
or gcd(P G (£ U A^)|q) = $ 7^ 5o which leads us back to case II. Now, it is known that 
if 5*0 is an initial solution of the Diophantine equation, S must have the general type: 
>S' = 5'o + Q*$, Qan arbitrary linear polynomial. These values of S ensure that $ is a 
common factor of P/|q and working as in the previous cases we can prove the last branch of 
the feedback-law design. In all of the above, causality is accomplished via a proper choice 
of the polynomial Q. 

Example 5.1 Lei us suppose that we have the nonlinear system: 

F = y{t) - 2y{t - 1) + 2y{t - l)y{t - 2) + 4y{t - 2)^+ ^y{t - 2)y{t - 3) - 2y{t - 2)y{t - 4)+ 

- Af + 2>u{t - 2) + l&u{t - if - lSu{t - l)u{t - 2) - 2u{t - lfu{t - 2)+ 
+5'ti(t - 2f + 2u{t - l)u{t - 2f - ]-u{t - 2f + u{t - lfu{t - 3) - u{t - l)u{t - 2)u{t - 3) + 
+ ^^t-2)V^-3) + 2/(t-lMi-2Mt-l) + ^y(t-2)Vi-l) + ^2/(t-l)|/(t-3Xt-l)- 
- 3fu{t - 1) - 4y(t - l)y{t - 2)u{t - 2) + y{t - 2fu{t - 2)+ 

+2y{t - l)y{t - 3)u{t - 2) - ^y{t - 3fu{t - 2) 

or using Se-polynomials F[y{t),u{t)] = 0, where F is the Se-polynomial A, appeared in the 
example 4-1- We want to find a feedback law u{t) = Sy{t) and some classes of desired sys- 
tems Ad, which are matched by the the closed-loop system,. To illustrate our approach we 
shall follow the MM- Algorithm in details. We shall follow the MM- Algorithm. In order to 
construct the sets Fi, F2, F3 we have first to calculate the sets of rules LV, GV, MV, DR, RV. 



means that any value of the parameters {wx^u,h,k, sx.cu.h.k) satisfies the condition ^ = 0. 
We use now the GCD- Values Algorithm to discover common factors of the members of 
C^. Let us further suppose that we are in step 1 of this algorithm and that n — 1, this 
suggests that we shall work with only one polynomial and all the others will be " elim- 
inated " by putting their companion coefficients equal to zero. Since ccoe/(Li^i q) = \ 
and ccoef{A) = 1, A = 6^ — 2di, this cannot be happened. Let us now suppose that 
n = 2 and we have again to work only with the polynomials -^1,1,0 and A by eliminationg 
the companion coefficients of all the others. A first glance at the coefficients indicates 
that the substitution wi, 1,0,3 = 0, wi, 1,0,2 = —4, 1^1,1,0,1 = 0, eliminates most of them, 
but finally the polynomial 1/6,1,0 will remain to the " game " with ccoe/(L6,i,o) = | md 
thus we have to go the n = 3 case. Here, we have a result. We obtain three polyno- 
mials with nonzero companion coefficients by means of the following rule of substitution: 
NZ = {wi,i,o,3 = 0,^1,1,0,2 = -4,Wi,i,o,i = 0,i(;6,i,o,i = 0,W6,i,o,o = 4}. These polynomials 
are A = {^4, Li,i,o, Le,i,o}. Now, we want to find a subset of NZ such that the polynomi- 
als A\nz = {^0 — 26i, —4^2 + ^4 + 'U'i,i,o,o^Oi 4(5o + ^2} have a common factor, but since 
gcd(5o ~ 26i,Ado + 61) = So this cannot be done and thus CF = 0. Working similarly, we 
explore all the different cases for various values of n and we finally conclude that always 
CF — and so GV — 0, which implies that Fi = 0. Let us now investigate the case F2. 
In order to get F2 7^ all the linear polynomials , including A and B, must have common 
factors, but since gcd(A, B) = Sq this is impossible and hence F2 = 0, too. R remains only 
the third case. This situation is similar with the previous one but without involving the no 
formal polynomials A and B. This will provide us with some degrees of freedom. Indeed, the 
application of the GCD-Values Algorithm will supply us with the following common factor 
for the polynomials of the sets C* and M. : ^ — —25q + 5i. Since Re,h — ^^^S h* ~ ^ 
finally get that $ G F3 7^ 0. It remains to specify the desired systems A^, which can be 
matched by the feedback-laws given in theorem 5.2. If A^ = then we seek for polynomials 
S and R such that 

6o-2Si = R* {-26o + di) + 3S2*S 

Solving this Diophantine equation we get a first solution Rq = — |5o + |5i and Sq = —^Sq. 
Therefore, the feedback-law upon request is u{t) = Sy{t), with S — —2So -\-di-\-Q* {—\So), 
Q arbitrary linear polynomial. This law matches the closed-loop system with Ad, for any 
set of initial conditions. If Ad has a common factor with $, then Ad must be a multiplier of 
$, since ^ is a prime linear polynomial we conclude that Ad = t*^, T a linear polynomial. 
In this case, the feedback-law u{t) = Sy{t), S as before, creates a closed-loop system which 
matches the desired system for any set of initial conditions {yo, yi), which satisfy the relation 
— 2t/i+yo — 0; (theorem 5.2). Finally, we choose Q properly so that causality to be satisfied. 
For instance, by taking Q — — 85o + 2^2; we get S — Si — or 

u{t)^y{t-l)-^y{t-2) 

Example 5.2 Let us consider a very simplified model for velocity control of an aircraft, 
[13], 

il^X2- f{xi), X2^-X2 + U, y^xi 

where Xi is the velocity, X2 is the engine thrust and f{xi) — x\ — 2xi is the aerodynamic 



continuous form: 



y + {2y- l)y + - 2y = u 



To discretize the above system we choose to use the Euler forward difference approxima- 
tions, i.e. y ^ y{w)-y{t) ^ - ^ -^[y{t) - 2y{t + 1) + y{t + 2)] and h ^ 1. A straight 
substitution and a proper shifting of the time delays, will give: 

—Vit) + y{t - 1) - - l)y{t - 2) + ^-y\t - 2) = -^-u{t - 2) (12) 

or using 5 and e-polynomials, F[y{t),u{t)\ =0 with F = \5q + 6i — "^8182 + +1^2 ■ We 
want to design a feedback law so that the closed-loop nonlinear system will match the stable 
linear system: 6y{t) = y{t — 1) + y{t — 2) or (— 65o + 5i + 82)y{t) — 0. We shall follow 
the MM-algorithm. The absence of nonlinear e-terms, as well as of 8e-terms, simplifies the 
complexity of the operations. Step 2 will give the set: C = {wio^o + i^'ii^i + 82, W2080 + 81, 
WS080 + 81, I ((-w^io + ^«n^30 - '^W2owIq - 2m;iim;2o^3o)^o+ {-2wioWu+2wIq + 2wiiwIq+ 
2uj1^uj3o - 4:W2oW2o - 4:Wu'W2oW3o)8i-\- {-2wio + 2^11^20 + 2w2o)82)}. At step 3 we form 
the following sets: LV — 0, GV — MV — {w^q — —2x,Wii — x — 2, W20 — —2,W3o — 
—2,x e R}, (we did not include all the operations in details for the sake of presentation) 
and we find a common factor $ = —28q + 81. Since gcd{Ad, $) = $ 7^ 80, theorem 5.2 will 
give S = Sq + Q * , where Q arbitrary, and Sq satisfies the equation R*^-\-S*B = A 
or R * {—28q + S * {—^82) = — + ^i- The latter equation will give S = — |5o 
and thus S — — |5o + Q * (— 25o + ^i)- -^'^ order to obtain a causal feedback connection we 
select Q = —^80 + 81 and finally S — — + 282. Hence, a feedback-law upon request is 
u{t) = — y — 1) + 2y{t — 2) . This will law will produce a close-loop system wich matches 
Ad for any set of initial conditions {yo, yi) which satisfy the relation —2yi +|/o- Please note 
that since the desired system is stable we have actually stabilized the nonlinear system, by 
using a feedback which is stable as well. 

6 Concluding Remark 

The aim of this paper was to describe an algebraic computational method for the study of 
a general class of non-linear discrete input - output systems that contain products between 
input and output signals. The model matching problem of these systems is solved through 
certain symbolic algorithms. The entire approach is based on a proper framework that 
involves the so-called (5e-operators and the star-product operation. We hope to be able to 
present current work on further applications of this method to other questions in a future 
paper. 
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